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ABSTRACT 

Our Galaxy is surrounded by a large family of dwarf galaxies of which the most massive are the 
Large and Small Magellanic Clouds (LMC and SMC). Recent evidence suggests that systems with 
the mass of the Local Group accrete galaxies in smaller groups rather than individually. If so, at 
least some of the Galaxy's dwarfs may have fallen in with the LMC and SMC, and were formed as 
part of the Magellanic system in the nearby universe. We use the latest measurements of the proper 
motions of the LMC and SMC and a multicomponent model of the Galactic potential to explore the 
evolution of these galaxy configurations under the assumption that the Magellanic system may once 
have contained a number of bound dwarf galaxies. We compare our results to the available kinematic 
data for the local dwarf galaxies, and examine whether this model can account for recently discovered 
stellar streams and the planar distribution of Milky Way satellites. We find that in situations where 
the LMC and SMC are bound to the Milky Way, the kinematics of Draco, Sculptor, Sextans, Ursa 
Minor and the Sagittarius Stream are consistent with having fallen in along with the Magellanic 
system. These dwarfs, if so associated, will likely have been close to the tidal radius of the LMC 
originally and are unlikely to have affected each other throughout the orbit. However there are clear 
cases, such as Carina and Leo I, that cannot be explained this way. 

Subject headings: galaxies: dwarf — galaxies: individual (Carina, Leo I) — Local Group — Magellanic 
Clouds 



1. INTRODUCTION 

Most of the mass in our Universe is believed to be 
an unknown, collisionless form of matter that can only 
interact gravitatiqnally with baryonic matter and itself 
(Efstathiou et al. 19901. The formation of large-scale 
structures in our Universe, in particular galaxy cluster- 
ing, is driven by gravitational forces exerted by this dark 
matter. Cold Dark Matter (CDM) cosmology holds that 
structures grow hierarchically, with small objects collaps- 
ing first and then merging to form more massive galaxies 
and clusters. The theory predicts that only dark matter 
halos with mass smaller th an approximately 10^ M© can 
form from 3tT fluctuations ( Read et al. 2006 ) in primor- 
dial density perturbations. Consequently, more massive 
systems can only form by subsequent accretion of these 
protogalactic fragments. For this reason, dwarf galaxies 
(dark matter subhalos in CDM cosmology) that contain 
luminous baryons and have not yet merged with a host 
galaxy, can, to some extent, be considered some of the 
most primitive buil ding blocks of our universe. However, 
it is now thought (Tosi 2003) that the inconsistencies 
between the observed properties of large galaxies and 



dwarfs are too many to believe that the former are built 
up only by means of successive accretions of the latter. 
We investigate two key issues identified with this model. 

The first issue has long been known, namely that the 
spatial distribution of the Milky Way satell ites shows 
asymmetric patterns and probably streains ( Kunkel & 
Demers||1976| |Lynden-Beil||1976[ |Hartwick'200(J) , partic- 
ularly when only the innermost satellites are taken into 
account (Figure [T]). The orbits of t hese satellites have 
been fou nd to be preferentially polar ( Zaritsky & Gonza- 



lez 1999 ) based on wealth of evidence including the a. 



fc Lynden- Bcll 1982 ), the orientation of the 
Stream (|Mathcwson et al 



ment of sate llites on the sky ( Kunkel fc Demers|1976 



distribution of sa .teliites ([ 
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1974 ), the three-dimensional 



revese et al. 1994) and their 



actual velocities ( Odenkirchen et ai.||1994 ). Knebe et al 



(2004) similarly find that nearby clusters such as Virgo 
and Coma possess galaxy distributions that tend to be 
aligned with the principal axis of the cluster itself. They 
conclude that either some dynamical process is responsi- 
ble or that the orbital parameters of the dwarf galaxies 
are imprinted on them at the time they enter the host 
halo. They concluded that this hypothesis can be ex- 
cluded at a 99.5% level, given the empirical constraints 
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Fig. 1. — Hammer-AitofT equal-area projection of the dwarf 
galaxy positions, with the 11 innermost Milky Way satellites high- 
lighted (black triangles). 



showing the Milk y Way potential to be spherical (Fell- 
hauer et al. 2006). Multiple arguments have been put 
torward to explain the disk-of-satellite s problem, primar- 
ily based on CDM simulation models. Kang et al. (2005 1 
argue that the observed distribution ot the Milky Way 
satellites is indeed consistent with being CDM subha- 
los, assuming that the satellites follow the distri bution 
of the dar k matt er within the Milky Way halo. |Zent- 
ner et al. (2005) used a semi-analytic model to iden- 
tiiy luminous satellites and showed that an isotropic dis- 
tribution is not the correct null-hypothesis; rather, the 
host halos are mildly triaxial (tending to be more prolate 
than oblate). Following similar lines to Libeskind et al. 
( [2005j , they found that the distribution of the Galac- 
tic satellite system is in fact consistent with being CDM 
substructure, albeit with a low probability. However, 
this disk-of-satellites is much less extended than nearby 
dwarf galaxy associations, and such an association falling 
into the Galaxy is unlikely to have produced the disk-of- 
satellites ( |Metz et al.||2009[ ). 

The second more specihc problem is that, if simple ac- 
cretion were the only mechanism for Galactic growth, 
then the spheroidal dwarf galaxy Leo I has an improba- 
bly high radial velocity. Mateo (1998) finds the systemic 
velocity of Leo I to be 287.0 ± 1.9 km s^^ for any sub- 
sample of the set of 33 red giants in the Leo I dwarf 
spheroidal galaxy for which they obtained radial veloci- 
ties from spectra taken using the HIRES echelle spectro- 
graph on the Keck telescope. This suggests that Milky 
Way satellites as a whole need to be considered in the 
wider context of the whole Local Group. Due to this ex- 
ceptional velocity, Leo I has an anomalously large effect 
compared to ot her satellites on estimat es of the mass of 
the Milky Way. Sakamoto et al. (20021 use Bayes' theo- 
rem and the assumption ot isotropic orbits for the satel- 
lites of our Galaxy to estimate a median probable mass 
of'--^2.5 X 10^^ M0. Interestingly, however, the exclusion 
of Leo I would lower this estimate to ~1.8 x 10^^ M©. 
Their result implies that Leo I may not be gravitationally 
bound to the Milky Way (in order to agree with other 
estimates of th e Milky Way mass derived by IW ilkinson] 
fc Evans| (|1999) and others, see Table jl]). 

In this paper we postulate an additional mechanism 
to explain the asymmetric dwarf galaxy distribution and 
the unusually high radial velocity of Leo I. We conjec- 
ture that some of the dwarf galaxies may have 'piggy- 
backed' in with the LMC-SMC binary pair, but had their 
bound orbits disrupted by interactions with the central 



Galactic potential. Dwarfs bound to the LMC initially 
will be closer than most dwarf galaxy associations, and 
the subsequent interactions may have resulted in a 3- 
body interaction that ejected Leo I and thus could re- 
produce its current high velocity while still explaining 
the narrow spacing of the disk-of-satellites association. 
To test this hypothesis we first explore the past orbits of 
the LMC-SMC under the effect of a central Galactic po- 
tential pair in accordance with the latest measurements 
of their proper motions. Numerically solving Newton's 
equations backwards in time, we proceed to find their po- 
sition in phase space at apogalacticon. We then model 
the infall of this configuration as an A^-body system con- 
taining the LMC and SMC as well as multiple dwarf 
galaxies (hereafter, all references to dwarfs and dwarf- 
dwarf interactions exclude the LMC and SMC) in bound 
orbits around the LMC. We evolve our system forwards 
in time and see whether it can reproduce the anomalies 
in the Local Group discussed above, and if so what was 
the state of the dwarfs when bound to the Magellanic 
system. 

2. FIDUCIAL MODEL 

The orbits of the LMC and SMC will have evolved over 
time as they interact gravitationally between each other 
and with the Galaxy over many gigayears. Assuming 
minimal perturbation of the orbits of the LMC and SMC; 
the orbits of dwarfs that accompanied the Magellanic 
system can be calculated from a previous apogalacticon 
to the present day and into the future. 

This is achieved by running a Monte Carlo suite of 
these multi-system models with varying initial conditions 
of the LMC and SMC (based upon the observed initial 
parameters, see [2.1). The LMC/SMC orbits are calcu- 
lated within a fixed Galactic potential and with gravi- 
tational interactions between the LMC and SMC (^2.2 1 



backwards in time for 3.5 Gyr and calculating the po 
sition of either the last or second last apogalacticon of 
the LMC. At this apogalacticon a number of d war fs are 
inserted in a bound orbit of the LMC (see ^2.3) and 
then traced forward in time to the present day and to 
+0.5 Gyr in the future, interacting with other dwarfs, 
the LMC/SMC system and the Galaxy. 

2.1. LMC and SMC orbital parameters 

The orbit of the LMC and SMC wiU be affected by their 
current positions, velocities and by their masses. The 
position and velocities used here for the LMC and SMC 
are given in Table [TJ Particular note should be taken of 
most recent meas urements of the proper m otions of the 
LMC and SMC by |Kallivayalil et al.| ( [20061 ), which give a 
relative velocity between the clouds at the current epoch 
of 105 ± 42 km s"^ These values imply that the LMC 
tangential velocity is ap proximately 100 km s^^ higher 
than previously thought ( Kallivayalil et ar]|2006 ). 

Given the low uncertainty m position, we assume that 
the LMC and SMC are at fixed locations given by the 
best estimates of observations. The velocity, with higher 
uncertainty, is randomly altered for each run assum- 
ing independent Gaussian uncertainties in the measure- 
ments, the line of sight velocity and the proper motions 
are calculated and then transformed to Cartesian veloci- 
ties around the Galaxy (for transform ations to Cartesian 
coordinates see van der Marel et~al][2002^ 
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TABLE 1 

Orbital parameters for the LMC and SMC adopted or derived from sources as referenced. 



LMC 



SMC 



References 



Line of sight velocity (km s ^) 
Proper Motions (W,N) (mas yr~^) 
Current Positions (a, S) (deg) 
Galactic Coordinates (1, b) 
Current positions (X, Y, Z) (kpc) 
Velocities (vx , vy , vz) (km s^-"-) 
GC radial velocities (km s^^) 
GC tangential velocities (km s^^) 



262.1 ± 3.4 
-2.03 ±0.08, 0.44 ±0.05 
81.9 ±0.3, -69.9 ±0.3 
280.5,-32.5 
-0.8,-41.5,-26.9 
3± 12, -268 ± 11, 252 ± 16 
89 ±4 
367 ± 18 



146 ±0.6 
-1.16±0.18,-1.17±0.18 
13.2 ±0.3, -72.5 ±0.3 
302.8, -44.6 
15.3,-36.9,-43.3 
-87 ± 48, -247 ± 42, 149 ± 37 
23 ± 7 
301 ± 52 



Harris &: Zaritsky 
ayalil ct al. f 





The masses of the LMC and SMC have large r uncer- 
tainties than that of the position and velocity, 'van der' 
Marel et al] (2002 ) obtained a mass for the LMC of 
(8.7 ± 4.3) X 10^ Mq within 8.9 kpc using an analysis 
of carbon stars. This mass is less than half that esti- 



mated by Schommer et al. ( 1992 1 who derive a mass of 
2 X 10^° M0 using radial velocities for several of the old- 
est star clusters in the LMC that lie well beyond 6 kpc 
of its center. We match the circular velocities of IvanI 



der Marel et al. ( 2002 1 with that of an Einasto halo at 
8.9 kpc to get a virial mass of 2 x 10^^ M^ for a halo 



virialized today. Using a simplified model of tidal radius 
^tid/^pcri = "itid/Mcaiaxy wc find that the LMC wiU 
have a mass of 5.6 x 10^° Mrr, within a tidal radius of 



30 kpc, only slig htly more massive than that of Schom- 
mer et al.Nl992L 



i-'or the SMU, [Hardy et al] ( |1989[ ) obtained a lower 
mass limit of 1.0 x 10^ Mi^ from observations of carbon 



stars, while Dopita et al. ( 1985 ) obtained 0.9 x 10^ from 
similar obs ervations ot planetary nebulae. More recently, 
Harris & Zaritsky (2006) determined a mass of 2.7-5.1 x 
10" M(7). We u se the velocity dispersion from Harris & 



Zaritsky (2006) to calculate an Einasto halo virial mass 
of 4.5 X 10^'^ M0, with a mass 5.7 x 10^ M© within a 
tidal radius of 19 kpc. 

These masses, somewhat above the calculated masses 
at smaller radii, also approximate the mass of a theo- 
rized co mmon halo b etween the LMC and SMC (see for 
example Bekki|2008 1 , without the problem of a large sub- 
halo witnin subhaio (in common halo models, the LMC 
halo will have a mass of ~20% of the common halo). 
These masses are also similar to those calculated from 
the st ellar mass content and seen inside A'^body simula- 
tions ( Boylan-Kolchin et al.||20TT ). 

For the model, we ignore any mass loss from the clouds 
(in particular that due to tidal stripping) , although mass 
loss is known to occur as evidenced by the amou nt of 



matter in the M agellanic Stream; (^2 x 10^ Mq Put 



man et ar]|2003 |. Dynamical friction between the LMC, 
SMC and the Calaxy is also calculated using the tidal 
masses of the LMC and SMC. The effect of dynamical 
friction between the LMC and SMC is ignored, although 
LMC/SMC binary systems are still foun d in -/V-body sim- 



ulations tha t include dynamical friction (Boylan-Kolchin 
et al.| 2011 ) suggesting that the survival time with dy- 



namical friction is sufficiently large for the Magellanic 
system to reach the present day with these masses. 

2.2. Potentials 



The LMC/SMC and dwarf galaxy orbits occur within 
the potential of the Galaxy and are altered due to inter- 
actions between the Magellanic clouds themselves and 
the dwarfs. 

There has been much work on Galax y potential models 
since t he fi rst basic schemes devised by Murai fc Fujimoto 
( 1980 1 and Fujimoto & Sofue ( 1976 1. In contrast to other 
recent approaches that take simple spherical Galactic po- 
tentials, we choose to model the Galaxy using a multi- 
component raodel of the Galactic potential developed by 
Flynn et al. (19961. This model matches known Galactic 
parameters such as the rotation curve, local disk den- 
sity and disk scale-length to high accuracy. The lack of 
noticeable warp in the Magellanic Str eam either in po si- 



tion on the sky or radial velocity — see Lin et al. ( 1995 ) 



suggests that the disk and bulge likely has small or neg- 
ligible effects on the orbits of the clouds, allowing us to 
have greater faith in our derived orbits. In this model the 
Galactic potential ^{R, z) is given in cylindrical coordi- 
nates, where R is the planar Galactocentric radius, and 
z is the distance above the plane of the disk. The total 
potential $ is modelled by the sum of the three different 
potentials: the dark halo ^h, a central component ^c, 
and a disk thus 



$7 



(1) 



The potential of the dark halo is assumed to be 
spherical and of the form 



1 



(2) 



where r is the Galactocentric radius (r^ = R'^ + z'^). The 
potential has a core radius tq, and Vh is the circular ve- 
locity at large r. The potential of the central component 
$(7 is modelled by two spherical components, represent- 
ing the bulge/stellar-halo and inner core components: 



GMci 



GAL; 



(3) 



where G is the gravitational constant, Md and rd are 
the mass and core radius of the bulge/stellar-halo term, 
and Mc2 and rd are the mass and core radius of the 
inner core. The disk potential is modelled using an 
analytical form t hat is a combination of thre e Miyamoto- 
Nagai potentials Miyamoto & Nagai (1975): 



D 



Dl 



D2 



D3 



(4) 
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TABLE 2 

Adopted parameters for the Galactic 
potential. 



Component 


Parameter 


Value 


Dark Halo 


ro 


8.5 kpc 




Vh 


220 km s-l 


Bulge/Stellar- halo 


rci 


2.7 kpc 




Mc, 


3.0 X 10^ Mq 


Central Component 


rc2 


0.42 kpc 




MC2 


1.6 X 10^° M0 


Disk 


b 


0.3 kpc 




Mb, 


6.6 X 10^° M0 




ai 


5.81 kpc 






-2.9 X 10^° Mq 




012 


17.43 kpc 






3.3 X 10^ M0 




0.3 


34.86 kpc 



TABLE 3 

Comparison of masses between various Milky Way models 
at radii of 50 kpc and 100 kpc. 



Model 



Af5o{10" Mo) A/ioo(10i2 M©) 



Adopted 

Wilkinson fc Evans| l |1999^ 
Kochanck (199(^1 
El & White ( 2 U08{ 
Smith et aL||j2UUY| 
SakamotoTt al. | | ["2002[ l 



4.3 

5.4 
4.9 



5.4 



1.4 

1.9 

1.81 
0.7-2.0 
1.8 



where 



1,2,3... (5) 



and the parameter b is related to the disk scale-height, 
a„ to the disk scale-leng th, and Mji^ t he ma sses of the 



three disk components. Flynn et al. (1996) discuss in 



detail the justification of this model and its close fit to 
observational data. We assume hereafter that it reason- 
ably accurately models the Galactic potential with the 
adopted parameters listed in Table [2] 

Solving Poisson's equation with the Galactic potential 
and integrating to 100 kpc gives a total mass of Mmw = 
1 .43 X 10 Mq . This value agrees well with Galactic mass 
derived in previous work for radii greater than 20 kpc (see 
Table |3|. 

Although the potential well of the Galaxy is by far the 
deepest, the LMC, SMC and dwarfs will interact with 
each other. We assume that the LMC and SMC (and 
each dwarf) consist of an Einasto halo, with a sharp 
boundary at the tidal radius. The acceleration experi- 
enced by one of these objects at xi by another at X2, 
separated by a distance r = ||xi — X2II is then ( |NicholS| 
fc Bland-Hawthorn|[2009l ) 



{3w2r,2-2-3/a exp(2/a)a~i+3/" 
X 7(3/a,2[r/r,]"/a)^^^i^ if r < rtidai,2 
'^^^r^-' ifr>nMai,2, 

(6) 

where Vs is the halo scale velocity, Ts the halo scale radius 



and a the Einasto scale factor (set here at a = 0.18). 

2.3. Dwarf Galaxy Orbital Parameters 

The relatively narrow spread in dwarf galaxy positions 
suggests that if dwarf galaxies did fall in with the Mag- 
ellanic system they would be more tightly bound than 
most dwarf associations seen today. We hence assume 
that the dwarfs were in a bound orbit of the most massive 
object in the Magellanic system, the LMC, at a previous 
apogalacticon of the LMC. 

Each dwarf is randomly assigned a periapsis and cir- 
cularity (and hence eccent ricity) according to the dis- 
tributions of Wetzel] ( [2011 ) extended down to the lower 
mass of the LMC. I hese distributions lead to a CDF for 
periapsis of 



F 



kpc 



:48.7{7-^(3/17,7„ 



' pen 



/r-vir.LMcl^-^^'bmax " 7min])}'/''-''(7) 

where, 7"^ is the inverse lower incomplete gamma func- 
tion and 



Tmin / 



where 



mill /max 



mm/ max 



7(3/17, [0.32r^i„/^,,/r™,LMc]"-''), (8) 
is the minimum and maximum allowed 



periapsis of the dwarfs. We choose a rmax of the tidal 
radius, as any dwarf that finds itself outside the tidal 
radius at closest approach will experience a greater at- 
traction to the Milky Way than the LMC and will hence 
be definitely lost. A minimum periapsis is chosen to be 
three times the exponential s cale length of the LMC dis k 
(disk scale length ^ 1.4 kpc van der Ma rel et al.| |2002), 
orbits close to the disk are likely to produce warping that 
is not observed within the LMC disk. 
The circularity CDF is given by 

F{ri) = 2.25r]^ "^ 2i^i(2.05, -0.66; 3.05; ?/), (9) 

where rj is the circularity of the orbit and 2^1(0, b] c; z) 
is the Gaussian hypogeometric function. 

Any orbits which would extend past the calculated 
tidal radius are excluded and a new one calculated. 

The position of the dwarf within the orbit is deter- 
mined by approximating the rosette orbit that arises 
from the spherical potential as a Keplerian orbit (the 
distributions also approximate it as Keplerian, but an 
extended NEW halo sees no significant change outside 
of low periapsis orbits, see Wetzel|2011 ) and is found by 
solving 



N 1 / 

P (r < a;rpo,.i) = - ( 2 



13 — arctan 



l-x(l-e)' 



HP) 



where e = y/l — rj'^ is the eccentricity of the orbit. 

By randomly determining the angular position of the 
orbit and the radial direction of travel, the position and 
velocity can then be determined in Cartesian coordi- 
nates. 

The mass of the dwarfs was set at either 1.1 x lO"^ M© 
within a tidal radius of 300 pc, consistent with the co m- 
mon mass of dwarf halos from Strigari et al. ( 2008 ) or 
a mass of 2.36 x 10* M© within the same tidal radius 
which represents a dwarf halo ten times more massive at 
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TABLE 4 
Group Infall Models. 







Mdwarf 


Initial 


Number 


% Unbound 


Model 




(lO'^ Mq) 


Apo 


of runs 


Today 


+500 Myr 


A 


1 


1.1 


1 


27800 


28% 


37% 


B 


1 


1.1 


2 


5600 


87% 


87% 


C 


4 


1.1 


1 


25200 


28% 


37% 


D 


4 


1.1 


2 


26300 


87% 


87% 


E 


4 


23.6 


1 


3900 


26% 


36% 


F 


4 


1.1 




2500 


8% 


9% 



Note. — N-q is the number of dwarfs originally bound to the LMC, 
Initial Apo is the number of apogalacticons before the present from 
which the orbits are tracked to the present day and the percentage 
unbound is the fraction of dwarfs that have enough energy to escape 
from an isolated LMC. Model F had no apogalacticon and was begun 
2500 Myr ago. 

the virial radius. The tidal radius of the dwarfs was set 
arbitrarily to enclo se th is mass, although any dependence 
on mass is small (j |3.1[ ). 

3. MODEL PARAMETERS AND THEIR EFFECTS 

In addition to the randomized position and velocity of 
the LMC/SMC and the dwarfs; the number of dwarfs 
bound to the LMC, the mass of these dwarfs, and when 
the dwarfs were bound to the LMC were altered. These 
parameters were altered throughout five different models, 
with the parameters shown in Table Hj Model F, which 
had no apogalacticon was calculated by multiplying the 
calculated error (see ^2.1| by a factor of three, in order 
to achieve the required velocities to be unbound to the 
Galaxy under this potential. 



3.1. The effects of the mass and number of dwarfs 

The amount of dwarfs will determine the amount 
of dwarf-dwarf interactions that occur, and within the 
comparatively small halo of the LMC these dwarfs 
could gravitationally disturb each other. The fre- 
quency of interactions between small dwarfs (with a mass 
A^dwarf /-^LMC ~ 10~^) bound to the LMC wiU be a fac- 
tor of their orbital distribution — if they themselves are in 
a disk-of-satellitcs around the LMC, not modelled here, 
interactions will be more frequent — and the amount of 
dwarfs present. 

The frequency and effect of these dwarf-dwarf interac- 
tions are examined by comparing the final energies of the 
dwarfs around the LMC in models A and C and models 
B and D. We consider a dwarf to be unbound at a point 
in time if it has enough energy such that it would es- 
cape the LMC considering only a two-body interaction. 
A two-sample Kolmogorov-Smirnov test could not rule 
out the null hypothesis that the distributions in final en- 
ergies from the LMC of models A and C or models B and 
D are consistent today (with probabilities that the dis- 
tributions are consistent of pxc > 0.3, pbb > 0.4) or at 
+500 Myr (pAC > 0.3, pbd > 0.9). This null resuh sug- 
gests that dwarf-dwarf interactions are rare with only a 
few dwarfs bound to the LMC, or that such interactions 
cause only slight perturbations in other bound dwarfs 
orbits. 

The effect of increasing dwarf tidal mass by a factor 
of ^20 is measured by the comparison of the energy of 
the dwarfs in models C and E. Applying a two-sample 
Kolmogorov-Smirnov test to the energy of the dwarfs 



from models C and E suggests the null hypothesis cannot 
be ruled out with probabilities of the same distribution 
being pcE 0.03 today and pcE > 0.2 at -fSOO Myr. 
The increase in p value at -1-500 Myr for all model com- 
parisons may be a consequence of the strongest LMC- 
Galaxy-dwarf interactions occurring near perigalacticon, 
approximately today, with any variance along the LMC- 
SMC orbit having a bigger effect now than in -f 500 Myr. 

3.2. Initial time effects 

Due to the large differences in the percentage of dwarf 
galaxies that become unbound between models that be- 
gan at the last or second last apogalacticon or that had 
no apogalacticon, see Table |4j it would be expected that 
there are significant variations in the energy distribu- 
tions of the dwarfs. Using a two sample Kolmogorov- 
Smirnov test shows significant differences between the 
final energy of dwarfs in models A and B today with the 
probability of the null hypothesis, that is the the final 
energies of dwarfs in models A and B coming from the 
same distribution, being pab < 10~^ with a KS statis- 
tic oi D = 0.70 and between models C and D today 
PCD < 10"^ with a KS statistic oi D = 0.69. These 
energy distributions are also significant between models 
C and F (pcF < 10"^, D = 0.20) and models D and F 
(pdf < 10^"^, D = 0.81). Similarly the differences are 
significant at -t-500 Myr {p < 10"^ between A&B, C&D, 
C&F and D&F). Even within bound dwarfs, the distri- 
butions in energy are significantly different {p < 10^^ A 
& B, C & D, C & F and D & F). As a run lost a higher 
percentage of dwarfs the remaining bound dwarfs were 
of lower energy; that is, bound dwarfs in model F were 
of higher energy than dwarfs in model C which were of 
higher energy than model D. This loss of high energy 
dwarfs, occurs as the LMC/SMC approaches perigalac- 
ticon. Runs in model F, which had orbits most altered 
by dwarfs due to the large timescales, were more likely 
to lose dwarfs the closer the perigalacticon passage (oc- 
curring approxiamtely today). 

4. RESULTS 

As the ranges of mass and number of dwarfs explored 
show little effect on the number of ejected systems, we 
focus on the orbits of the LMC/SMC and the orbits of 
the dwarfs in models C, D and F in order to examine 
these effects. 

4.1. Orbits of the LMC and SMC 

Due to the varying initial positions and velocities of the 
LMC and SMC, the orbits differ substantially between 
runs, with times for apogalacticons varying by hundreds 
of millions of years. Within model C, the SMC is bound 
to the LMC at the last apogalacticon in 13840 (55%) 
runs, within model D the SMC is bound to the LMC 
two apogalacticons ago in 2691 (10%) runs and within 
model F the SMC is bound to the LMC 2500 Myr ago in 
559 (22%) runs. A bound orbit for model D is shown in 
Figure [2] from t ^ —2500 Myr, with the current velocity 
for boththe LMC and SMC shown as a vector emanating 
from their current position. An unbound orbit for model 

D is shown in Figure [3] from t 3000 Myr, with the 

initial velocity again as a vector. Due to the large number 
of orbits in which the SMC may be unbound, this orbit 
is not necessarily typical of all runs. 
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Fig. 2. — An orbit where the SMC is bound to the LMC from model D with current position as circles and velocities shown as vectors 
(size of vector in kpc is the velocity in lO's of km s~^). The orbit of the LMC is shown as solid line and the orbit of the SMC is shown as 
a dashed line. 



If the SMC was unbound at the last apogalacticon 
while the dwarfs were bound, the dwarfs were more likely 
to be ejected {p < 10^'^), possibly as a result of the 
SMC and LMC becoming bound. However, this effect 
is small, with on average 0.08 more dwarfs ejected per 
run in model C when the SMC was unbound at the last 
apogalacticon, with no impact if the SMC was bound 
at the second to last apogalacticon (model D, p > 0.05). 
This suggests that the SMC becoming bound to the LMC 
may only remove dwarfs that would in any case become 
unbound at a future point. 



4.2. The orbits of dwarfs 

Despite being gravitationally bound at a previous 
apogalacticon, a large fraction of dwarfs become un- 
bound from the LMC, often ending up a great distance 
away from the Magellanic system at the present day. 

The final positions today (shown in Figure W] for model 
C, Figure [5] for model D and Figure |6] for model F) show 
clear differences due to the starting apogalacticon. 

For dwarfs that were bound at the last apogalacticon 
(model C), ~30% of these will be unbound today, how- 
ever, they remain relatively close to the LMC and are 
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Fig. 3. — An orbit where the SMC is not bound to the LMC from model D with current position as circles and velocities shown as vectors 
(size of vector in kpc is the velocity in lO's of km s"'^). The orbit of the LMC is shown as solid line and the orbit of the SMC is shown as 
a dashed line. 



spread out along the orbital path of the Magellanic sys- 
tem. These unbound dwarfs form a bimodal distribution 

leading and trailing the LMC. This bimodality illustrates 
the difBculty for a dwarf to remain unbound near the 
LMC, whose gravity will recapture most dwarfs close by. 

Of the dwarfs that were bound two apogalacticons ago 
(model D), a much larger fraction (~90%) have become 
unbound by the present day and have begun their own or- 
bits around the Galaxy. These unbound dwarfs, occupy 
three main types of orbits, a relatively small (^100 kpc) 
low eccentricity orbit around the Galaxy, a much larger 
(^300 kpc) highly eccentric orbit, with others remaining 



close to the LMC/SMC and travelling along similar or- 
bits. As in model C, dwarfs that are close to the LMC 

and unbound have a bimodal distribution leading and 
trailing, again showing that is difficult for a dwarf to 
remain unbound near the LMCs gravitational influence. 

Dwarfs that accompanied the LMC in an originally un- 
bound orbit (model F), very few have become unbound 
in the model today or at -1-500 Myr. These dwarfs are 
much more likely to have a significant effect on the or- 
bit of the LMC and SMC and correspondingly show a 
larger region of bound dwarfs as the model LMCs orbit 
is altered from the initial backwards integration. The 
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distribution of unbound dwarfs is similar, however, to 
those dwarfs bound at the LMC at the last apogalacti- 
con with a leading group and a slight overdensity trailing 
the LMC. 

Despite these dwarfs being ejected from the system, 
they remain in the same plane as the orbit of the LMC 
and SMC and within the disk-of-satellites seen around 
the Galaxy today, these unbound dwarfs also end up re- 
taining their original direction of travel despite poten- 
tially being ejected at high velocity in a complicated 
three-body interaction. All but one dwarf from models C 
and D are moving clockwise in the Y-Z plane, with the 
one dwarf travelling anti-clockwise having only a small 
velocity in the plane (uyz ^ 8 km s^^). 

Of note is that there are no dwarfs in models C, model 
D or F which comes close to the current position o f Leo I 
(-122 kpc,-117 kpc,189 kpc) ( [Kroupa et aD '2005). In 
addition with Leo I likely being unbound to the Galaxy, if 
it did infall bound to the Magellanic system an extremely 
rare multi-body reaction must have occurred to boost it 
away from the Galaxy. Simil ar ejections have b een seen 



is bound initially, see ^4.1) 



in large A^-body simulations (jSales et al 
satellites of the host system 



2007| ) between 
Leo i could therefore still 



have been associated with the Magellanic system, but 
unbound to the LMC when entering the Galactic halo. 

The large fractions of simulated dwarfs that were 
bound in the past and unbound today means that the 
lack of dwarfs observed to be orbiting the LMC today 
does not imply there were no dwarfs bound to the Mag- 
ellanic system in the past. 

4.3. Differences between bound and unbound dwarfs 

With a large fraction of dwarfs becoming unbound, 
many of which leave the Magellanic system, the prop- 
erties that allow dwarfs to become unbound can be ex- 
amined. In particular, the contribution of the orbits of 
the LMC and SMC and the initial orbit of the dwarfs 
themselves around the LMC. 

The effect of variations in the Magellanic system's or- 
bit on dwarf loss is examined by looking at how many 
dwarfs were lost in each run. If the LMC/SMC or- 
bits have significantly different loss rates, then future 
refinement of the orbit of the clouds may give evidence 
against dwarfs falling in with the clouds. If no physical 
property of the LMC/SMC orbits correlates significantly 
with dwarf ejection probabilities, then we can model ejec- 
tions as a Poissonian distribution. Hence if ~30% of all 
dwarfs were lost, as in model C, and with no contribu- 
tion from the orbits of the LMC and SMC, the number 
of runs which lose to 4 dwarfs would be a Poissonian 
distribution with ratios of (6748, 10530, 6162, 1602, 156). 
In fact the runs from model C lose these dwarfs in 
the proportions of (6866,10379,6137,1637,181), a dif- 
ference that is statistically insignificant (a x^-test gives 
p > 0.05). The effect of the orbits is also not seen af- 
ter the second to last apogalacticon in model D; again 
modelling as a Poisson process the number of dwarfs 
lost per run would be expected to be in the propor- 
tions (6,181,1895,8820,15397), while the results from 
model D give (6, 225, 1891, 8697, 15481); (a x^-test gives 
p > 0.01). This suggests that the positions of the dwarfs 
at apogalacticon are much more important than the or- 
bits of the LMC and SMC (excluding whether the SMC 



The impact of different orbits of dwarfs around the 
LMC at apogalacticon is also examined. Dwarfs which 
originally possess a low energy — that is, are deep within 
the potential well of the LMC — are less likely to be 
stripped either when the SMC passes close by or at 
perigalacticon, when the gravitational attraction of the 
Galaxy is largest. A Kolmogorov-Smirnov test compar- 
ing the initial energy of the dwarf orbits that remain 
bound and those that are unbound indicates that the 
distributions of initial energy are different for model C 
today and at -1-500 Myr (the null-hypothesis that the dis- 
tributions are the same has probabilities of ptoday < 10"'^ 
and p+500 Myr < 10^'^) and for model D today and at 
+500 Myr (ptoday < 10-^ p+500 Myr < lO'^). The frac- 
tion of dwarfs that remain bound as a function of initial 
energy is shown in Figure [T] Dwarfs that were bound 
one apogalacticon ago (model C) are still undergoing re- 
moval from the Magellanic system as it passes close to the 
Galactic center. Dwarfs that were bound two apogalac- 
ticons ago (model D) are no longer undergoing removal, 
with dwarfs that are bound now remaining bound over 
the next 500 Myr, with a handful being recaptured as 
they pass close by the LMC during their own orbit of 
the Galaxy. No matter when the dwarfs were originally 
bound, dwarfs that were deep inside the potential well 
originally are much more likely to be bound today and 
remain bound over the next 500 Myr, while dwarfs of 
energy above about —1.8 x 10^^ erg are all equally likely 
to be removed. 

The large tail of dwarfs with highly eccentric orbits 
seen in model D is also reflected in the energy-angular 
momentum diagram (Figure [s]) , with these dwarfs oc- 
cupying high angular momentum, high energy positions 
on this diagram. Relatively few unbound dwarfs occupy 
the same position on this diagram as dwarfs that end up 
bound, and any detected dwarfs that occupy this posi- 
tion will be less likely to have fallen in with the LMC on 
a previous orbit. 

5. DISCUSSION 

The appearance of a disk-of-satellites around the Milky 
Way, has been associate d with the infall of a group of 
dwarf galaxies (Li & White 200 8]) and in particular th e 
infall of the Magellanic system (D'Onghia fc Lake|2008 ). 
This scenario — an association of dwarfs falling onto the 
Milky Way — has been found to require a much more com- 
pact distrib ution than is seen in dwarf associations in the 
local group ( Metz et al.|2009 ). These bound associations 
were also unstable when orbiting a larger Galaxy. Few 
dwarfs who would spend their entire undisturbed orbit 
within the tidal radius remain bound after two pericentre 
passages, and almost half are lost within 500 Myr of the 
first pericentre passage. Here we examine the properties 
of dwarfs that were bound to the LMC in the past, and 
whether this scenario could be responsible for some of 
the dwarf galaxies seen today. As no dwarfs seen today 
are bound to the LMC, possibly excluding the SMC, any 
dwarfs associated in the past with the Magellanic sys- 
tem were likely to have been high in the potential well. 
If there were multiple small dwarfs bound to the LMC, 
they were unlikely to have interacted, with no difference 
seen between models which had only one dwarf associ- 
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FlG. 4. — The distance from the LMC of the dwarfs in model C at the present time and density contours of the final Cartesian location. 
The contours of dwarfs that are unbound are shown as (blue) solid lines and dwarfs which remain bound as (red) dashed lines. The contours 
give the fraction of all dwarfs within a 10 kpc radius in seven logarithmic steps from 10~^ to 10^^ with the contours increasing in line 
width with increasing density. Leo I is plotted for reference, but is off the plot in Z (at ~200 kpc). 



ated and those with four. The presence of smaller dwarfs, 
which possessed insufficient gas to form dwarf galaxies 
will be more numerous and may have an effect on the 
orbit of these dwarfs. This lack of interaction may arise 
from the LMC stabilizing the orbits of dwarfs against 
minor perturbations. Interactions with larger galaxies, 
such as dwarfs interacting with LMC sized galaxies, or 
the LMC/SMC interactions can produce noticeable dif- 
ferences in the orbits and morphology o f dwarf galaxies 
p'Onghia etaLl[2009] [Besla et al.|[201p| . 
(Jf the dwarfs listed by IKroupa et al. ( 



( 2005 1 as possi- 



bly part of the Magellanic system when falling in, only 
Sagittarius, Ursa Minor, Sculptor, Carina and Draco are 
located within the area that our model dwarfs end up. 
Sextans and Fornax, located slightly outside the range 
of dwarfs, may also have fallen in with the Magellanic 
system having become unbound before the p revious two 
apogalacticons. Using the velocities from |Lux et al.| 
( j2010| ) we examine the direction of travel of the dwarfs 
in the Y-Z plane. Ursa Minor and Sagittarius are trav- 
elling clockwise in the Y-Z plane consistent with them 
falling in with the Magellanic system, while the remain- 
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Fig. 5. — The distance from the LMC of the dwarfs in model D at the present time and density contours of the final Cartesian location. 
The contours of dwarfs that are unbound are shown as (blue) solid lines and dwarfs which remain bound as (red) dashed lines. The contours 
show the fraction of all dwarfs within a 10 kpc radius in seven logarithmic steps from 10~* to 10~^ with the contours increasing in line 
width with increasing density. The two unbound types of orbits for dwarfs are plotted schematically in the Y-Z plane. Leo I is plotted for 
reference. 



ing dwarfs have errors large enough to be consistent with 
both travelling clockwise and counter-clockwise. Of the 
remaining dwarfs, Carina is unlikel y to have intera cted 
with the LMC in the last 3 Gyr ( |Pasetto et al.|2011[ ) and 
correspondingly did not fall in with the Magellanic sys- 
tem. Future observations, with precise proper motions 
(e.g. from GAIA) may exclude more dwarfs using this 
simple rotation test, or via their positions on an energy- 
angular momentum diagram, with the dwarfs from our 
models ending up in only a few regions of the diagram. 



some of which are exclusively bound to the LMC (of the 
dwarfs seen today, only the SMC may be bound). 

The assumption of a spherically smooth and static po- 
tential for the Milky Way halo is known to be unreal- 
istic. A triaxial potential can have notic eable impacts 
on dwarf orbits over time ( Lux et aL]|2010 ), but this po- 
tential makes modelling the history ot observed dwarfs 
backwards much harder. On the other hand it may al- 
low more of the dwarfs observed today the possibility of 
being associated with the Magellanic system at an earlier 
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FlG. 6. — The distance from the LMC of the dwarfs in model D at the present time and density contours of the final Cartesian location. 
The contours of dwarfs that are unbound are shown as (blue) solid lines and dwarfs which remain bound as (red) dashed lines. The contours 
show the fraction of all dwarfs within a 10 kpc radius in six logarithmic steps from lO"'''^ to 10~^ with the contours increasing in line 
width with increasing density. The two unbound types of orbits for dwarfs are plotted schematically in the Y-Z plane. Leo I is plotted for 
reference. 



time. 

Models in which the Magellanic system is on its first 
orbit of the Galaxy show a much reduced loss. This 
is partially due to the perturbation of the LMC and 
SMC's orbits by dwarfs over long times. The possiblity 
of dwarf companions therefore needs to be taken into ac- 
count when determing the orbit of the LMC and SMC. 
The LMC, and to a lesser extent the SMC, were assumed 
to be tidally stripped before entering the Galactic halo. 
Tidal stripping that happens as the LMC enters the halo 



will increase dynamical friction and subsequently result 
in a lower velocity today than those present in model F. 
As most dwarfs are lost when the SMC binds to the LMC 
(an event that is likely to h ave happened after the LMC 



enters the Galactic halo; see Boylan-Kolchin et al. 2011 1 



or at perigalaction (approximately today), dwarfs in this 
scenario may be better represented by model C. 

The restriction that dwarfs initially spend their entire 
orbit within the tidal radius of the LMC (with respect to 
the Galaxy) clearly leaves out bound orbits (those with 
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Fig. 7. — The fraction of dwarfs that remain bound as a function 
of their beginning energy. A dwarf deep inside the potential well 
is unlikely to become unbound, as the orbit increases in energy up 
to ~ — 1.8 X 10-'-'' erg the dwarf is increasingly likely to become un- 
bound. Dwarfs that were bound two apogalacticoiis ago (model D) 
are just as likely to be stripped in +500 Myr as now, indicating that 
any stripping of the dwarfs from the Magellanic system will have 
already taken place. Dwarfs that were bound one apogalacticon 
ago (model C) are still undergoing stripping from the Magellanic 
system, with dwarfs more likely to be unbound at +500 Myr then 
the present day. 

negative energy with respect to the LMC). Dwarfs on 

these extra-tidal orbits are extremely likely to be tidally 
stripped and and could also be the source of some dwarf 
galaxies that reside in the disk-of-satellites. 



6. CONCLUSION 

We have modelled the orbits of a group of dwarf galax- 
ies bound to the LMC at a previous apogalacticon using 
a Monte Carlo approach, varying the orbits within ob- 
servational errors. Dwarf galaxies bound to the LMC 
are unlikely to interact with each other, and a significant 
fraction become unbound from the Magellanic system 
after only half an orbit. These dwarfs would be located 
around the Magellanic system today and likely be no- 
ticeably associated with it. Dwarfs that were botmd one 
and a half orbits ago — that is at the LMC's second-last 
apogalacticon — are over six times more likely to become 
unbound than remain bound, with many dwarfs being 
located either in an extended structure of orbits, or in a 
tight disk around the Galaxy. This disk encompasses the 
locations of a number of dwarfs observed today around 
the Milky Way, so a number of these dwarfs may origi- 
nally have fallen in with the Magellanic system and been 
captured by the Galaxy. The common rotation direc- 
tion of the dwarfs in this ring provides a test to rule out 
any counter-rotating dwarfs as originally associated with 
the Magellanic system. The extended disk-of-satellites 
cannot be explained by the dwarfs being bound to the 
LMC within the last two apogalacticons, and may have 
another origin. In addition, the anomalous velocity and 
position of Leo I is not explained by this mechanism, 
with no dwarfs in any of the simulations approaching 
the position, let alone the velocity, of Leo I. 
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Fig. 8. — Energy-angular momentum Diagram of the dwarfs in model D. Even though unbound and bound dwarfs are likely to be found 
in the same regions, they occupy different regions of the energy-angular momentum diagram. Unbound dwarfs which have highly eccentric 
orbits around the Galaxy have higher energy for the same angular momentum than dwarfs on low eccentricity orbits. The positions of 
the LMC and SMC on the energy-angular momentum diagram are plotted as a black square and circle respectively. The contours in the 
energy-angular momentum diagram show the fraction of all dwarfs within a 5000 kpc km s~^ X 5 ■ 10^^ erg g~^ box, in seven logarithmic 
steps from 10"'* to 10"^ with the contours increasing in line width with increasing density. The contours on the right show the fraction of 
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